Decoherence within a single atom 



Gyula Bene^ and Szabolcs Borsanyi^ 
^Institute for Theoretical Physics, Eotvos University, Pdzmdny Peter setdny 1/A, H-llll Budapest, Hungary 
^Department of Atomic Physics, Eotvos University, Pdzmdny Peter setdny 1/A, H-1117 Budapest, Hungary 

(February 9, 2008) 

An "almost diagonal" reduced density matrix (in coordinate representation) is usually a result 
of environment induced decoherence and is considered the sign of classical behavior. We point 
out that the proton of a ground state hydrogen atom can possess such a density matrix. This 
example demonstrates that the "almost diagonal" structure may derive from an interaction with 
a low number of degrees of freedom which play the role of the environment. We also show that 
decoherence effects in our example can only be observed if the interaction with the measuring device 
is significantly faster than the interaction with the environment (the electron). In the opposite case, 
when the interaction with the environment is significant during the measurement process, coherence 
is maintained. Finally, we propose a neutron scattering experiment on cold He atoms to observe 
decoherence which shows up as an additional positive contribution to the differential scattering cross 
section. This contribution is inversely proportional to the bombarding energy. 
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Decoherence is a remarkable quantum phenomenon 
which plays a central role in the (not yet fully under- 
stood) emergence of classical properties. It is also a prac- 
tically important effect in quantum optics, atomic physics 
and mesoscopic systems. Therefore, it can be of interest 
to discuss some of its less known aspects. In the present 
paper we show on a simple example that contrary to the 
common belief, 

1. an "almost diagonal" (in coordinate representa- 
tion) density matrix may be the result of an in- 
teraction with only a few degrees of freedom 

2. an "almost diagonal" density matrix does not nec- 
essarily imply the loss of coherence 

3. interaction may maintain rather than destroy co- 
herence. 

We also discuss the conditions when and how decoherence 
can be observed. 

Our primary example is a single ground state hydrogen 
atom. Later we shall also consider the case of a helium 
atom. 

The wave function of a freely moving hydrogen atom 
in the ground state may be written in coordinate repre- 
sentation in the product form 
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where 



R 
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stands for the coordinates of the center of mass and 
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for the relative coordinates. The vectors fp, fe refer to 
the position of the proton and that of the electron, respec- 
tively, while mp and me are the corresponding masses. In 
Eq.(0) 



Voir) 



1 



: exp 



a-B J 



(4) 



is the ground state of the hydrogen atom, where as = 
— ^ = 0.529 X 10~®cm stands for the Bohr radius. As 
the Schrodinger equation of the hydrogen atom separates 
in the R, r coordinates, the time evolution of the center 
of mass wave function does not influence the state ipo{r), 
and the product form (||) persists. Nevertheless, when 
expressing the total wave function in the original coordi- 
nates Tp, Te, the state of the hydrogen atom is not of a 
product form, instead, it is an entangled state. Indeed, 
both and foir) depend on both fp and fg. For sim- 
plicity suppose that AP as/?! <C mp/rrie (AP standing 
for the momentum uncertainty of the center of mass), or, 
equivalently, that 



AV < « We W 2 X 10^ — 
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{AV being the velocity uncertainty of the center of mass 
and Ve standing for the velocity of the electron), then 
ijjiR) ~ ip{'rp), and we arrive at the Born-Oppenheimer 
approximation |0] 



ipirp)ipo{re - rp) 
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Entanglement implies that the proton is not in a pure 
state, i.e., its density matrix p{fp, fp') is not a single pro- 
jector. Indeed, using expression (^ we get 
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Pirp,rp')= I (Pfe'i>irp,reW{rp',fe) (7) 
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where s = \rp — rp'l/as- Provided that the width of 
the center of mass wave function is much larger than the 
Bohr radius, the reduced density matrix (H) is "almost 
diagonal", i.e., its width (which is approximately 2 as) 
is much smaller than its length along the diagonal. The 
condition for this is AP as/h ^ 1, or 



Ay < 
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It is a much stronger condition than that of the Born- 
Oppenheimer approximation (cf. Eq.(^)). Let us empha- 
size that the " almost diagonal" structure in this example 
is due to the interaction with a single electron. 

An "almost diagonal" density matrix (in coordinate 
representation) is usually a result of "environment in- 
duced decoherence" and is considered the sign of classical 
behavior. In the present paper we consider the validity 
of this expectation and discuss the physical meaning of 
density matrices like 

It is of worth mentioning that we do not need to rely 
upon the Born-Oppenheimer approximation in order to 
arrive at an "almost diagonal" density matrix. That 
approximation is only convenient because it allows one 
to perform the integration in Eq.(^ explicitly, indepen- 
dently of the actual form of the center of mass wave func- 
tion. More generally, Eqs.(|l|), (Q) imply 
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This inequality shows that the offdiagonal matrix ele- 
ments of the reduced density matrix decay at least expo- 
nentially when increasing s = \rp — rp'\l as- 
ks an example, suppose that 
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where M 



7^e sa nip. Eq. ( pi] ) describes a 
Gaussian wave packet moving in free space. If condi- 
tion (||) does not hold, the integral in Eq.(^ cannot 
be expressed in a closed analytical form. Assuming 
the validity of (^, however, we obtain Eq.(^). We get 
a narrow density matrix (an "almost diagonal" one) if 



The deviation of the reduced density matrix of the pro- 
ton from the pure state can be characterized by Trp^. It 



is unity for a pure state and less than unity (but positive) 
for a mixed state. Using Eqs.([ll]), (||) we have 
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It is easy to see that in case of an "almost diagonal" 
density matrix (i.e., when z <C 1) Trp^ is much smaller 
than one (in fact, it goes to zero with z as z^), while 
in the other extreme of very large z it approaches unity, 
which corresponds to a pure state. 

As it is well known, a reduced density matrix gives 
all the possible information about the results of mea- 
surements done on the corresponding subsystem. Let 
us emphasize that such measurements are assumed to 
be infinitely fast, so that they give a snapshot about 
the state of the subsystem. It is obvious that coordi- 
nate measurements do not reveal the structure of the 
density matrix, as the probability distribution of the co- 
ordinates are given by the diagonal elements Pp{rp,fp) 
which is independent of the decay of the offdiagonal el- 
ements. Note that in the Born-Oppenheimer approxi- 
mation (cf. Eq.(^) Pp{rp,fp) = \ip{rp)\'^, which is com- 
pletely independent of the interaction between the proton 
and the electron. In contrast, momentum measurements 
are sensitive for the behavior of the offdiagonal elements 
of the reduced density matrix in coordinate respresen- 
tation. The probability distribution of the momentum 
is determined by the diagonal elements of the reduced 
density matrix in momentum representation, i.e., by 
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Introducing coordinates parallel and perpendicular to the 
diagonal of the reduced density matrix (in coordinate 
representation), i.e.. 
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we may write 

Pp{Pp,Pp) = 



ni = oirp + rp') 



rj_ = Tp - Tp 



^|dV^exp(-Jpp.n 



(15) 
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This last expression shows that the momentum distribu- 
tion is proportional to the Fourier transform of the ofFdi- 
agonal elements of the density matrix in coordinate rep- 
resentation, after having averaged them along the diago- 
nal. Obviously, the narrower the density matrix (in coor- 
dinate representation) becomes, the broader the momen- 
tum distribution will be. Using the Born-Oppenheimer 
approximation (|^) and the explicit expression (jl^), we 
get 
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where q — \pp ^ Polas/^ and zq = as/S. Note that 
the momentum distribution is independent of the time, 
which is a natural consequence of the fact that momen- 
tum (unlike coordinates) is now a conserved quantity. 

It is instructive to evaluate expression ( |i8| ) in two ex- 
treme situations: for 3> 5 and <C 5. The former 
case corresponds to a pure state of the proton, while in 
the latter case one has an " almost diagonal" mixed state. 
For a_B ^ ^ we get 
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It coincides with the momentum distribution of the cen- 
ter of mass, and is independent of the electronic motion. 
In the opposite case (which is our case of interest), for 
as ^ S we have 
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This coincides with the momentum distribution of the 
atomic electron (if q is identified with the momentum of 
the electron times as/h) and is independent of the center 
of mass motion. Let us emphasize that it is the momen- 
tum distribution of the proton, and the above coincidence 
is a consequence of momentum conservation and that the 
variance of the total momentum is negligible. Obviously, 
the width of distribution ( po|) is w fi/as, much larger 
than h/5, the momentum uncertainty corresponding to 
the state (pi]). This broadening marks the "almost di- 
agonal" structure of p{fp, Vp'). Below we shall point out 
that this sign of decoherence also shows up in neutron 
scattering cross sections which are experimentally more 
accessible. 

The above discussion raises the obvious question how 
one can then observe diffraction with atomic beams 
why interference is not destroyed by decoherence. In a 
typical diffraction experiment done with atoms the condi- 
tion of "almost diagonal" structure of /o(rp, rp'), Eq.(^ is 
indeed fulfilled. The essential point of the explanation is 



that the density matrix which characterizes a subsystem 
(in our case the proton) refers to instantaneous (infinitely 
fast) measurements. This means that one can observe de- 
coherence effects in a measurement done on the proton if 
during the interaction between the proton and the mea- 
suring device the state of the electron does not change 
significantly. As an estimate, one may compare the mo- 
mentum actually gained by the electron from the proton 
during the measurement to the momentum change of the 
electron that would take place if the electron followed 
the proton instantaneously (as in the Born-Oppenheimer 
approximation). One can also tell that decoherence ef- 
fects become observable if the measurement is so fast that 
the Born-Oppenheimer approximation for the electronic 
motion fails. Later we shall express the condition in a 
concrete situation quantitatively. The condition of ob- 
serving decoherence is obviously not fulfilled in a diffrac- 
tion experiment: the time of the measurement here is 
the time of fiight between the grid and the screen. Dur- 
ing this time the interaction between the proton and the 
electron maintains the form (|^) of the state of the atom, 
i.e., the Born-Oppenheimer approximation is valid all the 
time. For definiteness, let us consider a two slit experi- 
ment where the distance between the slits is much larger 
than as but smaller than 5. Be a{R,t) and (3{R,t) out- 
going waves emerging from the first and the second slit, 
respectively. Then the center of mass wave function can 
be written as 



V'(i?, t) = aa{R, t) + bf3{R, t) , 
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where a and b are the probability amplitudes that the 
atom goes through the first and the second slit, respec- 
tively. Eq.(^ remains valid until the atom reaches the 
screen, thus Eq.(|l|) is also valid during the time of flight. 
Initially, at the instant of time when the atom has passed 
through the slits a{R, t) and f3{R, t) are still narrow, sep- 
arated wave packets, but later on they broaden and over- 
lap. At the (only approximately determined) time in- 
stant to when the atom hits the screen Eqs.(|])-(||) are 
still valid, thus we get for the probability distribution 
P{rp) of finding the proton near a given point Vp the 
expression 



P{rp) = \aa{rp,to) + b(3{fpM)[' 
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This includes interference terms, too. Thus we see that 
despite of the "almost diagonal" density matrix (which 
has had such form already initially) coherence is not de- 
stroyed. Paradoxically, this is due to the same interac- 
tion (Coulomb attraction between the electron and the 
proton) which is responsible for the "almost diagonal" 
structure of the density matrix. 

In order to shed more light on the situation, let us 
consider what we would get if the interaction between 
the electron and the proton were turned off when the 
atom had passed through the slits. This can be done 
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in principle even experimentally by making use of the 
fact that during fast electronic processes the nucleus re- 
mains "frozen" . In molecular physics this is called the 
Franck-Condon principle Q. If the atom is irradiated by 
a suitable ultraviolet laser beam, it can be ionized. This 
is a fast process which practically does not influence the 
position and motion of the proton, but eliminates the in- 
teraction between the proton and the electron. In order 
to see the consequences let us express the initial state 
<I>(rp,re,i = 0) (i.e., still before the ionization) in the 
Schmidt (or biorthogonal) representation B, i.e., as 



$(fp, fe, i = 0) = ^ CjX]{rp)i 
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w aa{fp,t = O)ipo{fe - rsi) 
+bl3{rp,t = O)(poire - rs2) ■ 
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Here r^i, fs2 stand for the slit positions, Xj(^p)"S (i-e., 
a{rp,t = 0) and l3{fp,t = 0)) are the eigenstates of the 
reduced density matrix of the proton, while ^j(re)-s (i.e., 
'Poi'i^e — ^si) and (pQ{re " t^s2)) are the eigenstates of the 
reduced density matrix of the electron. The approximate 
orthogonality of (po(?\; — 'f^si) and foif,, — rs2) follows 
because |^^^ 1 — ^52! 3> as- During and after the ionization 
process the states of the electron-photon system evolve 
in time unitarily. The states Xji^^p) of the proton have a 
separate unitary time evolution. Therefore, we have 



\a\^\a{fp,t„)\^ + \b\^\/3{fp,t„)\^ , 
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i.e., interference terms are absent. It can also be un- 
derstood as a consequence of the fact that the "which 
way" information after the ionization can be obtained by 
measuring the electron without disturbing the proton. 

In the presence of the interaction, however, interfer- 
ence and coherence is restored. In that case, of course, 
separate unitary time evolutions for the proton and the 
electron do not exist, instead, they move together and 
the atom behaves (in free space or in slowly varying po- 
tentials) as a single unit. 

Let us return now to the question how one could ob- 
serve experimentally the consequences of the " almost di- 
agonal" density matrix (||). As mentioned above, deco- 
herence effects emerge if during the measurement done on 
the proton the interaction between the electron and the 
proton is negligible. One possibility, namely, photoion- 
ization behind the slit has already been mentioned. In 
this case the duration of the measurement is unchanged, 
but the electron- proton interaction is " turned off' . An- 
other possibility is to perform a very fast measurement 
on the proton. Below we suggest such an experiment. 

Let us consider low energy (a few eV-s) neutron scat- 
tering at a helium atom, prepared in a state where the 
width of the center of mass wave function is much larger 



than the atomic size. Low energy neutrons interact with 
the nucleus through a contact potential 



gS{fri - ra) 
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where g = 27r?i a/ ^, a standing for the scattering length 
and /i being the reduced mass of the neutron- nucleus sys- 
tem. There is also an interaction between the magnetic 
moments of the neutron and the electron, however, in 
case of the helium atom the contributions of the two 
electrons to the magnetic scattering cancel each other, 
provided that an inelastic scattering (excitation of the 
electrons to higher levels) is energetically not possible. 
This is why we need helium for this experiment rather 
than hydrogen. 

The condition for observing decoherence is that dur- 
ing the interaction time d/v {d standing for the nucleus 
size and v for the neutron velocity) the momentum trans- 
ferred to the electron 

d 1 2g2 
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is much less than the total change of the momentum of 
the electrons well after the collision, 2meV. (As before, 
nie and Ve stand for the electron mass and velocity, re- 
spectively.) Thus the condition is 



V > -1/ Ve « 4 X 10'' — 

V a_B s 
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In other terms, the bombarding neutron energy must be 
much larger than 0.08 eV . This ensures that the electron- 
a interaction is negligible during the neutron-a collision. 

The suggested measurement is just a measurement of 
the differential scattering cross section of the neutron- 
helium collision (the neutron is detected) under the con- 
ditions described above. We also assume that the in- 
coming neutron can be described by a plane wave in a 
satisfactory manner. First order time dependent per- 
turbation theory gives a suitable approximation for the 
wave function of the whole system, then we calculate 
the probability of observing an outcoming neutron at a 
given angle with arbitrary momentum. Thus we get for 
the differential scattering cross section the formula (in 
laboratory frame) 
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Here k, k' stand for the wave vectors of the incoming 
and the scattered neutrons, respectively. The explicit 
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appearence of the density matrix of the a particle is just 
a sign of decoherence, which in turn is a consequence of 
the neghgible interaction between the nucleus and the 
electrons during the collision. In order to calculate the 
density matrix of the nucleus of the helium atom we again 
apply the Born-Oppenheimer approximation (|^) where 
instead of the one-electron wave function we have a two- 
electron wave function. For the ground state we use the 
hidrogen-like wave function 



Vo(?'l,?'2,CTi,cr2) 
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na% \ as 

xx°(ai,a2). (28) 



Here Z* = ^ is the effective atomic number and 

16 

X (cri,cr2) stands for the singlet spin function. Using 
Eq.(p8|) we get for the density matrix of the nucleus 



Note that Eq.(p3|) is given in laboratory frame and the 
angular dependence of the leading order corresponds 
just to isotropic scattering in the center of mass frame. 
Decoherence shows up as an anomalous contribution to 
the differential scattering cross section that is inversely 
proportional to the bombarding energy. This contribu- 
tion is always positive and has maxima at = and {) = 
TT. At i<^„ = 1 eV these maxima are h{0)/q'^ ~ 8.2 x 10^'* 
and /i(7r)/g2 = 2.27 x IQ-^. 

In conclusion we add that other kinds of interactions 
(e.g., electron-photon or gravitational) may also result in 
an almost diagonal reduced density matrix. Again, one 
has to be careful in interpreting it as a loss of coher- 
ence. A more detailed discussion of these questions will 
be presented elsewhere. 
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Finally, by doing the integrations in the limit as ^ S 
and 



kaB — as > 1 



(32) 
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(note that Eq. (|26|) implies q ^ 3.4) we have (by taking 
rua = 4m„) 
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and 
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